The effect of electronic entropy on temperature peculiarities of the frequency 
characteristics of two interacting anharmonic vibrational modes in j3—Zr. 
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A 2D temperature-dependent effective potential is calculated for the interacting longitudinal and 
transverse L— phonons of f3 zirconium in the frozen-phonon model. The effective potentials obtained 
for different temperatures are used for the numerical solution of a set of stochastic differential 
equations with a thermostat of the white-noise type. Analysis of the spectral density of transverse 
vibrations allows one to determine the temperature at which /3-Zr becomes unstable with respect 
to the longitudinal L— vibrations. The obtained temperature value practically coincides with the 
experimental temperature of the f3 ^ a structural transition in zirconium. The role of electronic 
entropy in the /3— Zr stability is discussed. 

PACS numbers: 63.20Ry, OS.lOGg, 63.20 Kr, 71.15Nc 
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In recent years, structural transformations of zirco- 
nium induced by high temperatures and pressures have 
been intensively investigated both experimentally and 
theoretically. Considerable attention has been given to 
studying, at atmospheric pressure, the structural stabil- 
ity of the high- temperature /3— phase and the peculiari- 
ties of its transition into an a— phase existing at temper- 
atures below T — 1136K-. The commonly accepted view 
today is that the structural f3 ^ a phase transition at 
atmospheric pressure is closely related to the softening 
of the transverse phonon with wave vector k = ^[110] 
(iV— phonon) experimentally observed in /3— Zb^*^. The- 
oretical arguments in favour of this viewpoint are mainly 
based on calculations of the effective potential for the N- 
phonon performed in the frozen-phonon model^. From 
these calculations it follows that the effective potential 
of the N-phonon of Zr has a two-well shape, and, as a 
consequence, the phonon frequency squared in the har- 
monic approximation is negative. The imaginary phonon 
frequency points to the instability of /?— Zr in the ground 
state. 

Using the perturbation theory formalism for anhar- 
monic effects in crystals, it was shown^that at high tem- 
peratures the /?-phase of Zr becomes stable to the atomic 
displacements corresponding to the N-phonon due to its 
interaction with other transverse vibrational modes lying 
along the (110) direction. Using a modified pseudohar- 
monic approximationSi, we in turn succeeded in showing 
that in this case stability may be attained by merely al- 
lowing for the inherent anharmonicity of the N-mode&. 

However, the stability of /3— Zr as a whole is not de- 
termined only by the lattice stability to the atomic dis- 
placements corresponding to the N-phonon. The calcu- 
lations performed in the "frozen "-phonon modeli show 
that /3— Zr is also unstable with respect to the atomic dis- 
placements corresponding to the longitudinal vibrations 
with wave vector k — |(1, 1, 1) (L(-phonon). 

At pressures larger than the triple point (T — 973 K 
and P = 5.5 GPa) such displacements cause the tem- 
perature phase transition from the to the uj- phase 
(the A/_B2-type structure). At room temperature the 



cj-phase of Zr is experimentally observed in the pres- 
sure range from 2.2 to 30-35 GPa^'-. At zero pressure 
the instability of /?— Zr to this mode manifests itself in 
a sharp decrease of the longitudinal vibration frequency 
obtained from the spectra of inelastic neutron scattering 
near k — 1, 1)2.. To our knowledge none of the cur- 
rently available theoretical studies discusses the mech- 
anism of stabilization of the high-temperature /?— phase 
with respect to the Li displacements. Moreover, it is 
still unclear whether stabilization is due to the inherent 
anharmonicity of the L/ mode or to the phonon-phonon 
interaction. 

In Ref. we have suggested a model for two inter- 
acting modes with wave vector k — |(1,1,1) (a strongly 
anharmonic longitudinal L; mode, and an almost har- 
monic transverse Lt mode) embedded in a thermostat 
modelling their connection with the rest of the crystal. 

The temperature dependence of the vibration fre- 
quency of the Li and Lt modes was found by solving a 
system of stochastic differential Langevin equations with 
a thermostat of the white-noise type, which describe the 
motion of a particle in a two-dimensional effective po- 
tential W{x,y). The potential was calculated within the 
electron density functional theory for a series of simul- 
taneous atomic displacements x,y corresponding to the 
Li and Lt modes. The interaction of these modes was 
shown to result in an induced anharmonicity for quick 
transverse Lt vibrations which in the absence of interac- 
tion remain almost harmonic at all temperatures consid- 
ered. The presence of induced anharmonicity leads to the 
frequency broadening of the spectral density S (w) and to 
the appearance of a complicated temperature dependence 
of S{uj) of the Lt mode. 

The calculation of the probability density of the mean- 
square atomic displacements^fi shows that the probabil- 
ity of finding the atoms at the sites corresponding to the 
bcc lattice (/?— phase) increases with temperature. How- 
ever, even at temperatures much higher than the point 
of the /? — > a structural phase transition this probabil- 
ity is almost three times smaller than the probability of 
finding the atoms at the uj sites. In other words, even at 
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high temperatures the bcc lattice remains unstable with 
respect to the Li displacements. 

We have assumedifi the two-mode effective potential 
W{x,y) to be temperature independent. This is a rough 
assumption because, firstly, it does not allow for the 
crystal lattice expansion with increasing temperature. 
Secondly, at finite temperatures the potential W{x,y) 
should be calculated not from the total energy of a crys- 
tal in the ground state, E^i, but from the free energy 

F{T) = E,i~TSeU (1) 

where Sei is the electronic entropy and T is the tem- 
perature. In Ref. O the electronic entropy was shown 
to be of considerable importance in stabiHzing the high- 
temperature bcc phase of Zr. 

The effect of the crystal volume change on the two- 
mode potential W{x,y) will be detailed elsewhere. In 
this work we shall restrict ourselves to the calculation of 
the electronic entropy effect on W{x,y). We shall also 
find the temperature dependence of the spectral density 
of the longitudinal, L/, and transverse, Lj, modes with 
k = 2/3(111) and determine, at zero pressure, the stabil- 
ity region of the (3 phase of Zr with respect to the atomic 
displacements corresponding to these modes. To this end 
we shall successively consider the dependence of the elec- 
tron entropy Sei (T, x, y) and free energy F(T, x,y) on the 
temperature and atomic displacements x, y. The effective 
potentials Wt(x,?/) = F{T,x,y) ~ F{T, 0,0) calculated 
at each fixed temperature will be used in solving a set 
of stochastic differential Langevin equations of motion 
with a white-noise thermostat. An analogous stochastic 
approach was earlier usedi^ to model the lattice vibra- 
tions of a strongly anharmonic crystal. In Ref. 0, as 
distinct from this work, the dynamics of motion of a sin- 
gle longitudial Li mode in a one-dimensional temperature 
independent model potential W{x) have been studied. 

I. THE TWO-MODE EFFECTIVE POTENTIAL 

As in a previous papenifi, the unit cell of /3— Zr was 
chosen as a hexagonal lattice with three base atoms. 
A detailed description of the unit cell geometry, its re- 
lation with the bcc lattice parameters, and the atomic 
displacements corresponding to the vibrations with the 
chosen wave vector k — 2/3 [111] may be found, for 
example, in Ref. 0. The total energy was calcu- 
lated by the self-consistent full-potential LMTO (FP 
LMTO) methodiiii with the GGA approximation for 
the exchange-correlation potential termi^. The same set 
of FP LMTO parameters was used for all displacements. 
The one-centre expansions inside the MT spheres were 
confined to Imax = 6 . The MT sphere radii were chosen 
equal to Rmt — 2.20 a.u. Integration over k was per- 
formed on a (10,10,10) mesh equivalent to 166 fc-points 
in the irreducible part of the Brillouin zone. The total 
energy was calculated for 108 pairs of coordinates x, y 



(18 values for the variable x corresponding to the lon- 
gitudinal Li mode, and 6 values for y corresponding to 
the transverse Lt mode). The energies obtained were 
first approximated by a tenth-degree polynomial for the 
longitudinal vibrations and then by a fourth-degree poly- 
nomial for the transverse ones. 
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Lt - displacement (y) L| - displacement (x) 

FIG. 1: Two-mode effective potential for longitudinal and 
transverse atomic displacements corresponding to the vibra- 
tions with wave vector k = |(1, 1, 1) in /?— Zr. The displace- 
ments are given in units of the lattice parameter a. 

The two-mode effective potential calculated in the 
frozen-phonon model is shown in Fign](the potential part 
corresponding to the positive period of transverse vibra- 
tions is presented in the upper panel, while the lower pan- 
els show the potential cross-sections for the transverse 
{y — 0) and longitudinal {x — 0) vibration components.) 
As seen from the figure, the total energy is minimum 
at zero displacements (the atoms are locaHzed at the uj 
sites). Two local energy minima at a; w |0.166|,?/ = 
correspond to the bcc atomic arrangement. It should be 
noted that to obtain the oj structure observed in zirco- 
nium from the bcc lattice it does not suffice to merely 
displace the atoms along the x-direction, one should also 
increase the c/a ratio. Nevertheless, in this paper we 
shall use the term "uj phase" for a lattice with zero dis- 
placements. At y = 0, the effective potential for longi- 
tudinal vibrations has a three-well shape, which agrees 
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with the calculations in Ref. |j. Since the local minima 
are shallow, the bcc lattice of Zr is unstable with respect 
to small longitudinal vibrations. At large transverse dis- 
placements, the potential for the longitudinal component 
transforms into a two-well one. 

The left lower panel of FigQ] shows the cross-sections 
corresponding to the "pure" transverse atomic vibrations 
near the uj (solid line) and /3 (dotted line) centres. The ef- 
fective potential is seen to have in either case a parabolic 
shape, hence, the vibrations frequency for the transverse 
branch, in the absense of interaction with the longitudi- 
nal vibrations, must be described quite well by the har- 
monic approximation in both the uj and /? phases. 



II. CALCULATION OF THE ELECTRONIC 
ENTROPY 

The electronic excitation entropy term in (1) will be 
introduced as earlieuii. This means that our analysis 
is limited to the classical lattice-dymanics regime when 
the electron-phonon interaction effects become negligible 
and the electronic excitation entropy term Sei. is bare 
electronic entropy 



Se 



n[flnf - (1 - f)ln{l - f)]dE. (2) 



Here n{E) is the electron density of states (DOS), and 
f{E) is the Fermi distribution 



(3) 



(3 — l/(fc_BT). The chemical potential fJ.{T) is determined 
from the condition for the number of electrons z 



n{E)f{E)dE. 



(4) 



In paper ^3 Sei was calculated for the hep and hcc struc- 
tures of Zr and Ti. The electronic entropy was shown 
to be large and strongly nonlinear in temperature, as a 
result of both the volume and energy dependences of the 
DOS. 

In Fig. El n{E) is shown for different atomic displace- 
ments corresponding to the Li vibrations in the absence 
of transverse displacements (y = 0). 

The bottom picture corresponds to the w lattice (x — 
0.0), and the top one to the bcc lattice {x = 0.166). It can 
be seen that a considerable change in the electron spec- 
trum occurs at longitudinal atomic displacements. So, 
the DOS at the Fermi level, n{Ep), gets almost doubled 
on transition from the to to the bcc structure. The DOS 
at the Fermi level reaches its maximum at a; = 0.1, when 
the Fermi level passes through the peak of DOS. 

The electronic DOS at the Fermi level, n[Ep), and 
the entropy at T = BOOK calculated for the displace- 
ments corresponding to the L/ vibrations in /?— Zr are 
presented in FigEl It is seen that at room temperature 
the electronic entropy changes at Li displacements in the 
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FIG. 2: The DOS of Zr at different atomic displacements x 
corresponding to the L; vibrations. Zero displacements corre- 
spond to the Lj-phase. The energy is reckoned from the Fermi 
energy. 



same way as the DOS at the Fermi level. This result is 
in good agreement with the expression for the entropy of 
a strongly degenerate electron gas at low temperatures, 
known from the electron theory of metals^S 



-n[EF)T. 



(5) 



From this expression it follows that, at a constant tem- 
perature, the entropy is proportional to the DOS at the 
Fermi level. Nonetheless, it can be seen from the fig- 
ure that near x — Q and x = 0.16 this proportionality is 
violated. To clarify the situation, let us take a look at 
FigEI 

In Fig. Efa) the temperature dependence of the en- 
tropy is calculated for three atomic displacements corre- 
sponding to: (i) the uo lattice(a; = 0.0), (ii) the maximum 
DOS value x = 0.10, and (iii) x — 0.20. It is seen that, on 
the whole, in the temperature range considered the en- 
tropy is fairly well approximated by a linear function of 
temperature, which is in good agreement with the expres- 
sion for the entropy of a strongly degenerate free electron 
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FIG. 3: The DOS at the Fermi level (a) and entropy at T = 
300K (b) for the atomic displacements corresponding to the 
Li vibrations. 



gas J^. FigEfb) shows, however, that the coefficient of 
proportionaHty between entropy and temperature has an 
intricate temperature dependence. So, it increases with 
temperature at a; = 0.0, and diminishes with increasing 
temperature at a; = 0.10, while at a; = 0.20 the coeffi- 
cient grows with a rise in temperature up to T = 1200iir 
and decreases with further increase of T. In relation Q 
the coefficient of proportionality between entropy and 
temperature is merely the DOS at the Fermi level. The 
temperature dependence of S/T in Fig^l^b) can be ob- 
tained by replacated the Fermi energy in relation ^ for 
the chemical potential determined from the normalizing 
condition I^J- Indeed, as seen in FiglJl the DOS grows 
with increase in chemical potential for x = 0.0, and de- 
creases for X — 0.10. The intricate shape of the S/T 
curve at a; = 0.20 is due to the fact that with increasing 
temperature n{fj,(T)) passes through the maximum. Cor- 
respondingly, entropy grows at low temperatures and di- 
minishes at high ones. Thus, expression Q may be used 
to calculate the entropy in a wide temperature range. 

In Fig. the electronic entropy is calculated at T = 
1800K for the atomic displacements corresponding to 
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FIG. 4: The temperature dependence of entropy (a), and the 
coefficient of proportionality between entropy and tempera- 
ture (b) for different atomic displacements (a;=0, 0.1 and 0.2). 



both the longitudinal, L/, and transverse, Lt, vibrations 
with wave vector k — |(1, 1, 1) in /3— Zr. As seen from 
the figure, entropy is an intricate function of the dis- 
placements. The entropy minimum occurs at zero dis- 
placements corresponding to the atomic positions char- 
acteristic of the Lo structure (see the right lower panel). 

In the left lower panel of FigElthe changes in entropy 
are shown at transverse displacements for the uj (solid 
line) and bcc (dotted Hne) atomic configurations. The 
entropy is seen to increase for the atomic vibrations in 
the LU lattice and decrease in the bcc one. 



III. TEMPERATURE DEPENDENCE OF 
ENTROPY AND FREE ENERGY 

As seen from Fig. 31 the temperature dependence of the 
entropy differs for different atomic displacements x, y. In 
Fig. Elthe electronic entropy S(x, 0) is shown for various 
temperatures. The difference in entropy between the uj 
and bcc lattices is seen to increase with temperature, the 
entropy maximum shifting from a: = 0.10 to a; = 0.2. 

As the entropy enters in the free energy with neg- 
ative sign, such a change in entropy must lead to a de- 
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FIG. 5: Electronic entropy at T = 1800K for the longitu- 
dinal and transverse atomic displacements corresponding to 
the vibrations with wave vector k = |(1, 1, 1) in f3—Zv. The 
displacements are given in units of the lattice parameter a. 
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crease of the difference in free energy between the uj and 
P phases with increasing temperature, and, as a conse- 
quence, to the variation with temperature of the effec- 
tive potential shape. In the 2D case the free energy as 



a function of the displacements x, y has a shape similar 
to that in Fig^ Note that for each fixed temperature 
we obtain a particular effective potential WT(x,y). As 
an illustrative example, we show in Fig. [3 only the ID 
potential Wt{x, 0) for the L/ mode at different tempera- 
tures, as being the most important for further discussion. 
It can be seen that, at high temperatures, allowance for 
the electronic entropy results in a significant change of 
the effective potential. Although the potential remains 
strongly anharmonic at high temperatures, the energies 
of the bcc and uj lattices become comparable. 
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FIG. 7: Free energy for the longitudinal displacements at 
different temperatures. 



Thus, by calculating the free energy at different tem- 
peratures, we obtained the temperature-dependent effec- 
tive potetials for the interacting longitudinal and trans- 
verse vibrations of /? zirconium. In general, the magni- 
tude of the potential barrier between the oj and (3 phases 
decreases with increasing temperature, which should af- 
fect both the vibration frequency and the interaction be- 
tween the Ll and Lt modes. 



IV. EQUATIONS OF MOTION FOR A 
PARTICLE IN A POTENTIAL FIELD WITH 
RANDOM FORCES 

Let us consider the dynamics of motion of a pseudopar- 
ticle in the 2D effective potential W{x, y) shown in FigQ] 
As follows from the construction of W{x, y), the variables 
X, y are generalized collective variables corresponding to 
two waves of atomic displacements with the same wave 
vector A: = 1, 1) and different polarizations (longitu- 
dinal, X, and transverse, y). We assume that the effect of 
all the other phonons of the system may be represented 
as random forces acting on a system of two interacting 
oscillators. The dynamics of motion of two interacting vi- 
brational modes embedded in a thermostat with random 
forces is defined by the stochastic differentional equations 
of the Langevin type 
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cPx dW(x,y) dx ^ , . 
(Py dW{x,y) 



■ly- 



dt 
dy 



Fy{t). 



dy ' "dt 

Fx (t) , Fy {t) are random forces with correlators 

< F,{t) > = 
< F,{t)Fj{t') > ^ 2T-f5^j5{t-t'), 



(6) 



(7) 
(8) 



7 are coefficients of the vibration damping, T is the tem- 
perature of the thermostat. Assuming all the higher mo- 
menta to be zero aX n ^ 3 



The stochastic dynamical variables (coordinates X, Y 
and velocities Vx,Vy) found by solving Eqs.lO were 
used to calculate the autocorrelation velocity functions 
Ki{T) =« Vi{0)Vt{T) >>. To this end the whole inter- 
val of modelling time was divided into realizations. The 
values of dynamical variables X(to),Y(to) obtained at 
the end of each reaHzation were used as the starting val- 
ues for the next one. In the average, the length of one re- 
alization was Nr = 100000 steps (1.57 • 10-"sec). When 
calculating the correlators, the total number of reaHza- 
tions over which the averaging was performed amounted 
to 3000. 

The autocorrelation velocity functions are shown in 
FiglHlfor different temperatures. It can be seen that at 



< F{h)F{t2)...F{tn) >= 0, 



(9) 



we obtain the standard definition of the random process 
F{t) called Gaussian white noise. 

The coefficient in ^ was chosen so that Eqs.® de- 
scribe, in the limit of large times, the relaxation of the 
distribution function to the stationary Boltzmann distri- 
bution with temperature T 



P(X,v,i oo) 



exp 



mv'^/2 + W{X.) 



2hT 



(10) 



where X(t) is a dynamical vector variable, and the mean- 
square velocity v in the limit t ^ cxd is equal to the 
thermodynamically equilibrium value at temperature T 



= lim < V 

t — 'OO 



>= hT. 



(11) 



The equilibrium distribution lll()|l results from the ac- 
tion of two opposing tendencies. Owing to the pres- 
ence of random forces F the velocity V (and, hence, 
X) becomes a random value, while the term describing 
the damping 7V suppresses V, tending to reduce it to 
zero. The set of equations @ was solved by numeri- 
cal integration of stochastic differential equations, using 
a method^ being a generalization of the Runge-Kutta 
scheme to the stochastic differential equations. We used 
a 4- step method of third order with the parameters from 
Ref. 
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FIG. 8: Autocorrelation velocity functions < K(0)K(t)) > 
(upper diagrams) and < Vy(0)Vy(t)) > (lower diagrams) at 
different temperatures T. 



When solving the set of stochastic differential equa- 
tions, a particular effective potential W(T) was con- 
structed for each fixed temperature T in (HJ . The calcula- 
tion was performed with a time step At = 1.57-10~^^sec., 
the number of steps in the realization being Nsh = 3 ■ 10^. 
Thus, the total time of modelling was ~ 5 • 10~*sec., 
which is vastly greater than the period of vibrations of 
the chosen modes. The test calculations for Ng^ = 9 • 10^ 
have shown that the result remains uchanged with such 
an increase of the modelling time. 



low temperatures the autocorrelation functions oscillate 
for a time that exceeds the period of atomic vibrations, 
whereas at high temperatures the correlations decHne 
rapidly. This means that at high temperatures the sys- 
tem quickly forgets about its initial state and comes into 
a stationary one. The temperature dependence of the cor- 
relation time is determined by the conditions imposed on 
random external forces F{t) acting on the system of two 
oscillators (see relations II7I8|I ). 
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The first relation reflects the fact that the mean force 
acting on the system of oscillators is equal to zero. The 
second condition impHes that the interaction with the 
thermostat is practically instantaneous, and the suc- 
cessive interactions do not correlate with one another. 
As seen from Eq.JHl, the magnitude of random exter- 
nal forces increases with temperature, which results in a 
faster damping of the correlators. It follows from Fig|H| 
that at low (T = 300K) and high(T = 1800K) temper- 
atures the autocorrelation functions have practically the 
same oscillation frequency, while at intermediate temper- 
atures there are additional oscillations. This is especially 
noticeable for the autocorrelation of the y component of 
the velocity at T = 1200K, and for < Vxit)Vx{to) > at 
T = 900K. 

The spectral vibration density Si{uj) found from the 
autocorrelation velocity function Ki(r) 



Si{uj) 



cos(tJT)Ki(r) dr 



(12) 



is presented in Fig. ^a) for different temperatures. The 
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FIG. 9: Spectral density of the L— phonon calculated for the 
temperature dependent (a) and temperature independent (b) 
efTective potential. 



first maximum of the spectral density at frequencies close 



to w w ITHz is due to the longitudinal vibrations, and 
the second one, near oj « 4THz is connected to the trans- 
verse vibrations. As seen from the figure, the tempera- 
ture variation leads to a significant change in the spectral 
density of both longitudinal and transverse virations. 

We first consider the fine structure of the spectral den- 
sity of longitudinal vibrations, Sx{uj), and its evolution 
with temperature. At low temperatures there is only one 
peak of frequency ujf « l.QTHz, which corresponds to the 
atomic vibrations near the sites appropriate to the cu lat- 
tice (a: = 0). With increasing temperature, the frequency 
of such vibrations diminishes, which is due to a slight 
anharmonicity of the effective potential (FigO near the 
origin. At temperature T = 600K, overbarrier vibra- 
tions arise in the system, leading to the appearance of 
an additional peak in the vicinity of wf « 0.9THz in the 
Sx{(^) curve. With a rise in temperature, the probability 
of overbarrier vibrations increases, and at T = 900K the 
share of such vibrations becomes dominant. This results 
in an increase of the spectral density intensity Sx{cd) at 
Lu fa 1.0 THz and in its decrease at lu > 1.5 THz. This 
tendency persists with further increase in temperature, 
the maximum of the spectral density shifting to greater 
frequencies. As seen in FigEta), at T = 1800K practi- 
cally all Li vibrations are the overbarrier ones. 

An analogous change in the spectral density is also ob- 
served for the transverse atomic displacements (the fre- 
quency range w w 3 — 4.5 THz). Again, at room temper- 
ature we see one maximum of frequency ujf w 4.2 THz, 
corresponding to the vibrations near the sites of the uj 
lattice. Then, with increasing temperature, an additional 
peak appears in the range ~ 3.5 THz. With further 
rise in temperature this peak shifts to higher frequencies 
and becomes the main one. 
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FIG. 10: The temperature dependence of the frequency tJmaa; . 
The temperature-dependent effective potential is shown by 
the solid line, the temperature-independent one by the dashed 
line. 



In FigElthe soHd line shows the temperature depen- 
dence of the vibratonal frequency at which the spectral 
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densities Sx{uj) and Sy{uj) have a maximum value (indi- 
vidual values of ujmax were estimated for the longitudinal 
and transverse components of S{lu)). As seen, for both 
types of vibrations there is a discontinuity in ujrnax ■ From 
the above it follows that the discontinuity occurs at Ttr, 
at which the intensities are equal: Sx{u!f) = Sx{(^2) or 

For the longitudinal vibrations the temperature T^^ lies 
in the range from 600 to 900 K, and for the transverse 
ones Tj^ « 1200 K. The fact that the temperatures T/^ 
differ is due to different nature of the two-hump structure 
of Sx{ijj) and Sy{uj). In the former case, the presence of a 
structure with similar intensities of the spectral density 
reflects the fact that the overbarrier vibrations equal in 
number the vibrations near the lo centres, while in the 
latter case the similar intensities indicate that the shares 
of vibrations near the bcc and lj sites are equal. There- 
fore, the spectral density of the longitudinal vibrations, 
Sx{u}), carries no immediate information about the pro- 
portion of time the system spends respectively in the bcc 
and the uj lattice. Such an information is contained in 
the spectral density of the vibrations, Sy{oj). It is seen 
from the above that at a temperature close to 1200K the 
probabilities of finding the system in the uj and the bcc 
phase are equal. The temperature T^J. can be consid- 
ered as the temperature of equilibrium between the two 
phases with respect to the atomic displacements charac- 
teristic of the longitudinal L; vibrations with wave vector 
k = 2/3(111). 

For reference, in FigEtb) is plotted the spectral density 
of vibrations calculated for a temperature-independent ef- 
fective potential. In this case for all temperatures the 
effective potential was the same, equal to the poten- 
tial found from the total energy without allowance for 
entropy. A comparison between FigEl^a) and FigEtb) 
shows that the greatest discrepancy is observed for the 
transverse vibrations. Namely, the two-hump structure 
of the spectral density of transverse vibrations persists 
up to T = 18007^, which far exceeds the phase transi- 
tion temperature in zirconium. And at T = 1200K the 
main contribution to the spectral density comes from the 
vibrations with frequency lo w 4.2THz, corresponding 
to the vibrations near the equilibrium position of the uj 
structure. The calculated values oiojmax for this case are 
shown in Fig^jby a dashed line. As seen from the fig- 
ure, taking account of the temperature dependence of the 
effective potential leads to a decrease in the frequency of 
both longitudinal and transverse vibrations in the system 
and to a fall in temperature T^^. 

Thus, the calculation performed shows that the al- 
lowance made for the electronic entropy when modelling 
strongly anharmonic lattice vibrations results in a de- 
crease of the calculated temperature value at which the 
l3 phase of zirconium becomes stable to the atomic dis- 
placements corresponding to the L phonon. 



VI. CONCLUSION 

In this work we have used f3 — Zr to discuss the ef- 
fect of temperature on the effective potential W{x, y) and 
the spectral density of vibrations of the longitudinal, Li, 
and transverse, Lt, interacting modes with wave vector 
k — 2/3(111). The effective potential acting on these 
modes at zero temperature was calculated in the frozen- 
phonon approximation within the electron density func- 
tional theory by the FP LMTO methodii^. For nonzero 
temperatures, the potential Wt was defined as the dif- 
ference of the free energies Ft calculated for different 
atomic displacements, WT{x,y) = FT{x,y) — Ft(0, 0). 
The required electronic entropy was defined through the 
density of electron states and the Fermi-Dirac distribu- 
tion function. 

The results obtained show that in zirconium in the 
temperature range up to 2000K the entropy may be well 
approximated by a linear function of the density of elec- 
tronic states. The shape of the effective potential W{x,y) 
for the Li and Lt vibrations strongly depends on the 
temperature. With increasing temperature the height of 
the energy barrier between the bcc and uj structures de- 
creases. As a result, the bcc lattice of Zr becomes stable 
with respect to these modes at lower temperatures. An 
analysis of the vibrational frequency of the transverse Lt 
mode showed the temperature, at which the (3 phase of 
Zr becomes unstable with respect to the longitudinal L; 
displacements, to be Tl — 1150 ±25K. This value practi- 
cally coincides with the experimental temperature of the 
(3 ^ a transition: T/j^a = 1136K«i. Thus, the bcc lattice 
of Zr is dynamically stable with respect to the longitudi- 
nal Li vibrations owing to the intrinsic anharmonicity of 
this mode. Of considerable importance is the electronic 
entropy which substantially enlarges the stability region 
oi (3 — Zr with a decrease in temperature. It is commonly 
believed today that under atmospheric pressure and low 
temperatures the structural instability of /3-Zr is due to 
the softening of the transverse phonon with wave vector 
k = ^(110) (iVt-phonon)2*ifii. The result obtained in this 
paper points to the fact that at zero pressure there is 
at least one more process responsible for the structural 
instability of bcc zirconium: the atomic displacements 
corresponding to the longitudinal Li phonon. This brings 
up the question as to whether the coincidence of the tem- 
peratures T^^ and T^^q is accidental, or the L; mode is 
really of considerable, if not decisive, importance in the 
/3 — > a structural transformation in zirconium. 
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